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We derive a number of exact relations between equilibrium and nonequilibrium quantities 
for spin glasses in external fields using the Jarzynski equality and gauge symmetry. For 
randomly-distributed longitudinal fields, a lower bound is established for the work done 
on the system in nonequilibrium processes, and identities are proven to relate equilibrium 
and nonequilibrium quantities. In the case of uniform transverse fields, identities are proven 
between physical quantities and exponentiated work done to the system at different parts of 
the phase diagram with the context of quantum annealing in mind. Additional relations are 
given, which relate the exponentiated work in quantum and simulated (classical) annealing. 
It is also suggested that the Jarzynski equality may serve as a guide to develop a method to 
perform quantum annealing under non-adiabatic conditions. 
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1. Introduction 

The remarkable recent developments in nonequilibrium statistical physics, the Jarzynski 
equality^' and fluctuation theorem, ^"^^ have given a paradigm to understand dynamical 
behavior through a direct relationship between a nonequilibrium process and equilibrium 
states (Jarzynski equality) or the probabilities of a nonequilibrium process and its inverse 
process (fluctuation theorem). The purpose of the present paper is to examine dynamical 
properties of spin glasses by use of the Jarzynski equality and the fluctuation theorem and 
establish exact identities and inequalities involving equilibrium and nonequilibrium quantities 
for spin glasses. 

Properties of spin glasses have been studied extensively for many years by experimental, 
numerical, and analytical methods. Although most of the theoretical aspects have been 
understood fairly satisfactorily at the mean-field level, it is still difficult to establish ana- 
lytical results for finite-dimensional systems. One of the exceptional successful approaches is 
an analysis by the gauge symmetry.^*^' Not only equilibrium quantities but also nonequi- 
librium behavior can be evaluated without approximations.^^) Our previous study^^^ has re- 
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vealed several new types of exact equalities between nonequilibrium processes and equilibrium 
properties of spin glasses by Tise of the gauge symmetry in addition to the above-mentioned 
nonequilibrium relations. In that paper, we have considered an application of the method 
of gauge transformation to the Jarzynski equality to shed new light on the possibilities to 
use a relatively unconventional sampling method, annealed importance sampling or its im- 
provement known as the population annealing,^^"^^) in the theoretical studies of spin glasses. 
These methods use a relation analogous to the Jarzynski equality while changing the tem- 
perature similarly to simulated annealing^^) and show outstanding performance comparable 
to the exchange Monte Carlo method. However, our previous study treated the change in 
temperature and fell short of the investigation of the performed work corresponding to the 
direct change of the parameters in the Hamiltonian. In the present paper, we revisit the prob- 
lem of nonequilibrium processes in spin glasses and develop a theory to analyze the work done 
to the system by the change of the strengths of external fields. The resulting identities and 
inequalities represent simple relations between equilibrium and nonequilibrium quantities in 
spin glasses in external fields, which should be added to the list of rare exact results on static 
and dynamic properties of spin glasses 

This paper is organized as follows. In §2, we analyze the work performed on spin glasses 
by randomly-distributed longitudinal fields. We then discuss the possibility of non-adiabatic 
quantum annealing starting from equilibrium ensembles and analyze the work performed on 
spin glasses by the transverse field. In addition, we give highly non-trivial relations between 
two completely different processes, simulated and quantum annealing. The last section will 
summarize this paper. 

2. Spin glass in random fields with gauge invariance 



In the present section we discuss the ± J Ising model of spin glasses in random fields on 
an arbitrary lattice, 



where erf is the z component of Pauli matrix at site i and is considered in this section to be a 
classical variable taking the values ±1. The distribution function of quenched randomness of 
Jy is specified as 



2.1 Model 




(1) 



P{J,j) = p6{Jij - J) + (1 -pMJij + J) 



2 cosh /3pJ 



(2) 
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The parameter /3p has been defined by e^'^^P'^ = (1 —p)/p to recover the middle expression 
of eq. (2) from the right-most expression for Jij = J and Jij = —J. The function h(t) in eq. 
(1) represents the time-dependent strength of the field starting from h{Q) = ho and evolves 
toward h{T) = Ht, and r)i is the quenched random variable following the distribution function 

F(,.) =<,%.- 1) + (1 + 1) = 5^2^. (3) 

The parameter hp is given by exp(— 2/3p/ip) = (1 — g)/g. The following analyses can readily be 
applied to other distribution functions of Jij and hi as long as they satisfy a certain type of 
gauge symmetry. 

The above Hamiltonian is invariant under the gauge transformation given by the following 
simultaneous changes of the signs of interactions, random fields and spin direction: 

Jij Jijiiij, Vi VS, cfi cfiii (Vi,j)> (4) 

where is a gauge variable taking ±1. While the Hamiltonian (1) does not change under the 
above gauge transformation, the distribution functions (2) and (3) are modified as 

P(T ^- ^M^^p-hjUj) p. ^ _ exp{::i/,^,iuCi) 

^ 2cosh^pJ ' ^^''^ij- 2coshPphp ■ ^ ^ 

These properties help us to derive the results in the following sections in addition to the exact 
value of the internal energy and several exact/rigorous results on a subspace of the phase 
diagram known as the Nishimori line (NL), fi = Pp and h(t) = hp,^'^'^^^ where /3 is the inverse 
temperature. 

2.2 Jarzynski equality 

The Jarzynski equality is useful to relate equilibrium and nonequilibrium processes. The 
Jarzynski equality in the present problem for a specific configuration of {Jij} and {r/j} is 
written as 

/^-|3W^ Z^{hT;{Jij},{Vi}) ff.. 

Zf3{ho;{Jij},{m}) 

where {Jij}, {Vi}) is the partition function of the system (1) (when the system is kept 

in equilibrium with the Hamiltonian of eq. (1)). The symbol (• • • )rf denotes the average over 
all possible nonequilibrium processes that start from equilibrium with inverse temperature /3 
at t = and evolve following the master equation. The Jarzynski equality holds independently 
of the details of dynamics and the functional form of h{t). 

We use here a discrete time representation for simplicity, t = kAt (A; = 0, 1, 2, ■ ■ • ), where 
At is a small time unit. Correspondingly, the spin configuration changes as {crjo {t = 0), 
{cr}i {t = At), ■ ■ ■ , {a}n {t = nAt = T). Notice that each {cr}fc stands for an instantaneous 
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classical configuration of N spins as = {cr|^(t), cr|(t), • • • ,af^{t)} at time t = kAt. The 

performed work in the short interval At is given by the energy difference due directly to the 
change of the Hamiltonian 

-SWk = {h{{k + l)At) - h{kAt)} Vi<^l (7) 

i 

The total work W is given by the summation of this infinitesimal work over time, J2k ^^k- 



2. 3 Work performed by random fields 

The Jensen inequality reduces the Jarzynski equality (6) to the well-known expression of 
the second law of thermodynamics, 

(W^)rf > -\ log ( = nhr, {J.,}, {m}) - F{ho; {J.,}, {v^})) . (8) 

(3 Zf:j{ho;{Jij},{r]i}) \ ) 

The quantity in the large parentheses is the difference of the equilibrium free energies for the 
initial and final Hamiltonians. It is difficult to obtain the explicit form of the free energy for a 
specific configuration of {Jij} and \r]i\. Therefore we usually evaluate the configurationally- 
averaged quantity over all realizations of {Jy} and \r]i\. Let us thus consider the configura- 
tional average of both sides of the above inequality as 

where the square brackets with the subscript ^p, denote the configurational average fol- 
lowing the distribution functions of {Jij} and {r/j}, eqs. (2) and (3). As proved below, the 
right-hand side of eq. (9) reduces to, under the NL condition Pp = P and hp = ho, 

mM,M £ '^m^^VM) - fog (^) , (10) 

where N denotes the number of spins and -D(/3, /io|/9, hr) is the Kullback-Leibler divergence, 
D{p,ho\P,hT)= E mho; {J.,}, im}) log ^f'\\Yx irl! • (11) 

Here we defined the probability for the configuration of {Jij}, {iji} summed up over all possible 
gauge transformations, 

where A^^ expresses the number of bonds, and expressions in eq. (5) have been used. Since 
the Kullback-Leibler divergence is non-negative, the work performed by random fields during 
a nonequilibrium process from the NL condition (i.e. the left-hand side of eq. (10)) does not 



log 



Zpihr; {Jij}, {vi}) 



ZB(ho;{Ji4},{rii})_ 



(9) 
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> -- 



(13) 



To prove eq. (10), we apply the gauge transformation as 
Z0(hT;{Jij},{rii])' 



log 



Z^{hQ;{Jij},{qi}) 



= E 



exp {I3p JiMj + f^pf^P rii^i) Zp{hT] {Jij], {m}) 



(2 cosh PpJ)^B (2 cosh ^php)^ 



log 



Zp{hQ;{Jij],{r]i}) ■ 



(14) 



{Jij} Am} 

Then we sum it up over all possible configurations of {^j} and divide the result by 2^ to 
obtain 

Zp{hT;{Jij},{m})' 

Zi3{hT;{Jij},{rii}) 



log 



= E 



Z0{ho;{Jij},{Vi}) . 

Zpp{hp;{Jij,{Vi}}) 



log- 



(15) 



2^ {2 cosh PpJ)^B ^2 cosh Pphp)^ " Zf^iho; {Jij}, {%}) ' 

According to eq. (12), this leads to eq. (10) under the NL condition Pp = P and hp = ho. 

Instead of the inequality (8), we may consider the configurational average of eq. (6) itself. 
Then we obtain the following simple equation by similar calculations as in the above analysis, 

AT 



RF 



/ cosh/3/iT 
/3,ho \cosh/3/io 



(16) 



Again, by the Jensen inequality, a lower bound of the performed work is evaluated, leading 
to the same inequality as eq. (13), 



N 



f cosh/3/iT 



(17) 



This is a bound looser than the previous one, eq. (10), since the KuUback-Leibler divergence 
does not appear here. We emphasize that eq. (16) is highly non-trivial because it gives an 
explicit bound for the work performed on spin glasses in nonequilibrium processes. 



2.4 Relation with different quenched randomness 

By the method used in the previous study on the Jarzynski equality in spin glasses without 
external fields, we can derive several exact relations between equilibrium and nonequilibrium 
quantities for spin glasses in the presence of random longitudinal fields. The results in this 
section are relatively straightforward generalizations of those in our previous paper, and 
therefore detailed derivations are omitted here. 

It is not difficult to show an identity for an observable Ot measured at the final step of 
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Ore 



{0)p,hT 



Zi3{hT;{Jij},{rii}) 



(18) 



RF ' Zjs{hQ;{Jij},{'qi}) ' 

where {0)i3^hT denotes the equilibrium (thermal) average of the observable O in the final 
equilibrium state with the strength of the random field hx- This relation leads us to several 
unexpected equalities in spin glasses. Calculations by use of the gauge symmetry, similar to 
the above one, give the following identities 

/ cosh (3hT 



i^RF(T)e-^^) 



RF. 



m 



RF) 



cosh f3hQ 



N 



(19) 



We remark that diS'erent concentrations of quenched randomness, represented by ho and 
hx, appear in both sides on the equality. The same relation holds for other gauge invariant 
quantities. 

In addition to the gauge-invariant quantities, the gauge transformation yields exact equal- 
ities on the local magnetization cr? and correlation function aga^, both of which are gauge 
non-invariant, as 



(aar)e-^'^)RF 



(20) 



and 



(ao^(r)a,^(r)e-/3^)RF 



'',ho 



/ cosh Phr 
\ cosh l3ho 



N 



(21) 



These results relate physical quantities measured in quite different environments. They imply 
the possibility that equilibrium physical quantities in spin glasses, the right-hand sides of eqs. 
(20) and (21), can be evaluated from nonequilibrium calculations (left-hand sides) in different 
parts of the phase diagram with the aid of annealed importance sampling or population 
annealing. 

Furthermore, the fluctuation thcorcm^^^^ yields the following relation for a quantity that 
depends on the intermediate spin configurations, 0{{a}o, {cr}i, • • • , {ctIt), as 



0{{(T}o,{a}i,--- ,MT)e 



RF 



Z,3{hT:{Jij},{Vi}) 

Z(s{ho;{Jij},{vi}) ' 

(22) 



where Or denotes the observable that depends on the backward process hr ho, and the 

angular brackets with the superscript 'rev.' express the nonequilibrium average for the back- 
ward process. An application of the gauge transformation to eq. (22) gives the following exact 
equality for the autocorrelation function O{{a}o, {cji, • • • , {o'}t) = 



a!{0)a!{T)e 



I RF. 



[{<yM<yi{T)WWr ^^^^ 



N 



(23) 
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This result states that the equihbrium autocorrelation function in the reversed process can 
be computed from the nonequilibrium process in a different part of the phase diagram. 



3. Jarzynski equality for quantum annealing 

3.1 Spin glass in transverse field 

Let us next consider another system with a transverse field, 

i?QA(i) = -9{t) Jv^^^J - (1 - 9{t)) To (24) 

where g{t) = t/T, which changes from to 1 as t goes from to T. This system is used in 
quantum annealing for search of the ground state of the spin glass Hamiltonian^^"^^) 

Ho = -Y.Jv^t^h (25) 

{ij) 

which is shared with eq. (1). The whole Hamiltonian is invariant under the gauge transfor- 
mation, 

(7f ^ af , (7f ^ Cif^f, (^t ^ ii(^t. Jij ^ JiMj (^hj), (26) 

where ^i{= ±1) is a gauge variable. Notice that this transformation is designed to preserve 
the commutation relations between different components of Pauli matrix. 

The adiabatic theorem guarantees that a sufficiently slow decrease of the strength of the 
transverse field (i.e. large T) changes the trivial initial state, the ground state of — ToX^jCrf, 
to the nontrivial ground state of the target Hamiltonian Hq. This is a special case of quantum 
annealing, quantum adiabatic computation.^^) Quantum adiabatic computation, however, is 
known to be unable to solve efficiently certain instances of hard optimization problems. ^^'^^^ 
Thus, instead of the adiabatic control, we analyze a method to repeat non-adiabatic quantum 
annealing (small or intermediate T) starting from a state chosen from equilibrium ensemble, 
not necessary the ground state. We may not be able to easily reach the ground state of Hq 
by such processes even if we start from a very low-temperature state since the system does 
not trace the instantaneous ground state as in the adiabatic evolution. We instead need to 
repeat the process many times to accurately evaluate the average of the exponentiated work 
over non-adiabatic processes appearing in the Jarzynski equality. In this way, the problem of 
long annealing time is replaced by a problem of very many repetitions of non-adiabatic (pos- 
sibly quick) evolution. We analyze such non-adiabatic quantum annealing using the Jarzynski 
equality and gauge symmetry. 
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3.2 Non-adiabatic quantum annealing 

Initially we pick up a state from the canonical ensemble for Hq\{0) = —^oYli'^f ^i^d 
then let it evolve following the time-dependent Schrodinger equation. The performed work 
in the present quantum problem is given by the difference between the outputs of projec- 
tive measurements of the initial and final energies, W = Em{T) — En{0). Here m and n 
denote the indices of the instantaneous eigenstates measured at the final and initial steps, 
HQA{T)\m{T)) = EmiT)\m{T)) and i/QA(0)|n(0)) = En{0)\n{0)), respectively. The Jarzyn- 
ski equality is^^'^''') 

where Z{t; {Jij}) is the partition function for the instantaneous Hamiltonian (24). The left- 
hand side of eq. (27) expresses the average of the exponentiated work over all realizations of 
non-adiabatic processes starting from the equilibrium ensemble. 

Following the prescription of the Jarzynski equality, we consider a repetition of non- 
adiabatic quantum annealing starting from the equilibrium ensemble. The initial Hamiltonian 
is given only by the transverse field, which means a trivial initial distribution. Consequently, 
the Jarzynski equality (27) for non-adiabatic quantum annealing reduces to 



- (2cosh/3roW (28) 



3. 3 Work performed by the transverse field 



Let us take the configurational average of eq. (28) over all realizations of {Jij} for P = Pi 
and = as 

^ Zp,{T-{J,,}) 



(2 cosh /3i To) 



N 



The right-hand side is written explicitly as 



(e-'^^^)QA 



exp (/32 Efe) «^*^) %(T;{J,,}) 



(29) 



(30) 



f-^ (2 cosh ^2 J)^^ (2cosh^iroj 
Let us apply the gauge transformation and sum over all possible configurations of the gauge 
variables {^i}- We then obtain, after dividing the result by 2^, 



(e-^^^)QA 



/32 



fr^. 2^(2 cosh M^'' (2 cosh ^iFq)^ ' 



(31) 
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Fig. 1. Two processes of non-adiabatic quantum annealing in eq. (33). The left-hand side of eq. (33) 
corresponds to the annealing process ending at the upper-right black dot and the right-hand side 
terminates at the lower- left dot. Three phases (F: Ferromagnetic, P: paramagnetic, and SG: Spin 
Glass) are separated by solid curves and a vertical line. The dotted line expresses the NL. 



A similar average of the exponentiated work on spin glass with the inverse temperature /32 
and the parameter for the quenched randomness /3i gives 



/3i 



{■hj} 



Zp,{T-{J,,})Zp,{T-{J,,}) 
2^(2 cosh /3iJ)^s (2cosh/32ro)^' 



(32) 



Comparing eqs. (31) and (32), we find the following relation between different non-adiabatic 
processes, 

f cosh. j3iJ\^^ ( cosh/32ro^ ^ 



/32 



/3i \ cosh /32 J y Vcosh/SiFo 



(33) 



Figure 1 describes the two different paths of non-adiabatic quantum annealing related by 
this equality. If we set /32 = in eq. (33), (implying p = 1/2, the symmetric distribution 
or the high-temperature limit), we obtain a simple equality on the performed work during 
non-adiabatic quantum annealing 



(e-^i^) 



QA 



(cosh/3iJ)^B 



(34) 



" (coshAFo)^' 

The symmetric distribution (/32 = on the left-hand side) makes it possible to reduce the 
right-hand side to the above trivial expression. It is remarkable that non-adiabatic quantum 
annealing, which involves very complex dynamics, satisfies such a simple identity irrespective 
of the speed of annealing T. 
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Similarly to the classical case discussed in §2.4, we can formulate a relation between 
equilibrium and nonequilibrium quantities as 

(0.e-n«AHO),|£iMI, (35) 

where Ot is an observable measured at the final time. The angular brackets with subscript (3 
is the equilibrium (thermal) average by the final Hamiltonian H(^js^(T) = Hq. If we choose Hq 
as the observable O and take the configiuational average with the symmetric distribution for 
both sides of the above equality, we obtain 

/n (T\.-PW\ 1 JiVB(cosh^J)^ 

(^QA(r)e )qaJ ^ = — (eosh^ro)^ ^^^^ 

The quantity on the right-hand side is obtained with the relationship given by the gauge 
transformation as 



Zp{T,{J,,}) 
;^'^(2cosh/3ro)^ 



(cosh /3J) 



B 



„= (coshgro)" 

The identity [{Hq)^]^ = -JAT^tanh/^J on the NL^^'^^) leads us to eq. (36). Equation (36) 
reads 

[(^QA(T)e-/^^)QA]o ^ , „^ 

^ = -NbJ tanh^J. (38) 

[(e '''^)qaJo 

Similarly to the classical case, we have obtained the equilibrium quantity (right-hand side) 
defined away from the part of the phase diagram where the original non-adiabatic process is 
performed (left-hand side). 



3.4. Exact relations involving inverse statistics 



Let us next take the configurational average of the inverse of the Jarzynski equality, eq. 
(28), as 

r 1 n \ (n^^c\^ /Q-p_\-^ 

(39) 



1 




"(2cosh/3ro)^ 


L(e-^^)QAj 




_ZfiiT-{J,^}) 



When the gauge transformation is applied to the right-hand side, we have 



1 



QA 



E 



cxp %T.{ij) JiMj) (2cosh/3ro)^ 
(2cosh/3pJ)^B Zp{T-{,hj}) 



(40) 

/3p izcosnp,.--" ■ 

By summing the right-hand side over all possible configurations of {^j} and dividing the result 
by 2^, we find 



1 



L(e-^^)QA, 



{Jij} 



Zp^XT;{'h,}) (2cosh/3ro)^ 
2^(2 cosh PpJ)NB z^{T; {Jij}) 



(41) 
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Fig. 2. Two different nonequilibrium processes of non-adiabatic quantum annealing are related 
through eq. (43). The same symbols are used as in Fig. 1. The white circle denotes the target 
of the process on the right-hand side of eq. (43), whereas the black dot is for the left-hand side. 



Under the NL condition /?„ = /3, this equation reduces to 



1 



.(e" /QA 
Comparison of eqs. (34) and (42) reveals 



(cosh /SFf 
(cosh^J)^£ 



QA 



-I3W) 



QA 



(42) 



(43) 



As depicted in Fig. 2, two completely different processes are related by this equation: One 
toward the NL and the other for the symmetric distribution. 

Let us further consider the inverse of eq. (35) for Ot = ^^fc^j- We take the configurational 
average of both sides under the NL condition as 



(af^^e-/^^)QA 



(cosh (3Tq 
(cosh (3 J) 



Nb 



1 



(44) 



The quantity on the right-hand side is unity by the gauge transformation as has been shown 
in the literatures.^*^' We thus obtain a simple exact relation 



1 



«^|e-^^)QA 



N 



/3 



(cosh /SFq) 
(cosh/3J)^B' 



(45) 



which is another exact identity for processes of non-adiabatic quantum annealing. 



3.5 Possibility as a solver 

The last part of this section is devoted to a discussion on the possibility to measure 
equilibrium quantities through non-adiabatic quantum annealing. The ratio of eqs. (28) and 
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(35) gives 



(e-''^)QA 



(46) 



This equation suggests that the equihbrium (canonical) average under the Hamiltonian Hq on 
the right-hand side can be estimated by the non-adiabatic quantum anneahng on the left-hand 
side. This fact may be useful in the evaluation of equilibrium average when it takes a very 
long time to equilibrate the system in Monte Carlo simulations as in spin glasses, since the 
left-hand side is evaluated without slow adiabatic processes. Nevertheless we should be careful 
because the average on the left-hand side involves a non-extensive quantity, the exponentiated 
work, whose value fluctuates significanlty from process to process: Remember that the average 
on the left-hand side is calculated by many trials of annealing processes. Thus, rare events 
with large values of the exponentiated work (i.e. l3\W\ 1) would contribute to the average 
significantly, and we have to repeat the annealing process very many times in order to reach 
the correct value of the average. 

The low-temperature limit of the above argument gives us a suggestion for non-adiabatic 
quantum annealing to identify the ground state. If the temperature is lower than the energy 
gap between the ground state and the first excited state of Hq and /9 S> Tq, the Jarzynski 
equality (27) is reduced to 



where EQs{{Jij}) is the ground-state energy of Hq. This equation suggests that the ground- 
state energy can be evaluated by the repetition of non-adiabatic quantum annealing, which 
implies that we would hit the correct ground state if we repeat the non-adiabatic process many 
times. We have thus replaced the problem of long annealing time in quantum adiabatic com- 
putation by another problem of many repetitions of non-adiabatic (fast) quantum annealing. 
It usually takes very many, typically exponentially many, repetitions to correctly evaluate the 
left-hand side of eq. (47). Thus the difficulty has not been relaxed in general, but the present 
new perspective may lead to different methods and tools than conventional ones to attack the 
problem. 

A related remark is that, for classical systems, an improvement by a modification of the 
dynamics has succeeded in estimating the free energy difference through fast nonequilibrium 
processes through the Jarzynski equality.^^^ A similar idea has been realized for quantum 
annealing using another degree of freedom.^^) 



(e-^'^)QA « exp{-pEGs{{Jij}) + PNTo) 



(47) 
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It is possible to establish identities to relate completely different annealing processes. 
The Jarzynski equality holds for simulated annealing of the classical system Hq through the 
'pseudo work' defined by 

-SY{kAt) = - (^(fc+i)At - PkAt) Ho, (48) 

where the inverse temperature is changed from (5kAt to /3(^k+i)At at the A;th step. Then the 
Jarzynski equality is 

where Y{/3o — >■ /3t) is the sum ^^(^y(A;At). Using the pseudo work Y in the weight of 
importance sampling, population annealing is implemented to produce the equilibrium en- 
semble. 



4-1 Identities for simulated annealing 



Let us consider the configurational average of eq. (49) over quenched randomness of {Jij} 
in Hq with j3p = j3r, where j3r is a given constant, 

exp {^T.{ij)J^) Z[s^{T;{Jij}) 



(e 



-Y{I3o-^I3t) 



)SA 



E 



(2cosh^,J)^s Z^,{T;{Jij}) 



(50) 



We apply the gauge transformation and sum over all configurations of {^i}. Division by 2^ 
of the result leads to 

Zf3^{T-{J,,}) Z^^{T;{J,,}) 



SA 



/3r 



E 2^(2 cosh /3,J)^s Z0,{T;{Jij}) 



(51) 



It is useful to write a similar equation for the case starting from Pq and ending at f3r with the 
quenched randomness characterized by Pt, 

,-y(/9o-i-/9r)\„.l _ ^/3t(^; Z(^^{T;{Jij}) 



(e- 



(52) 



(53) 



f^. 2^ (2 cosh /3r J)^« Zf,, (T; {Jij}) ' 
Equations (51) and (52) relates two different annealing processes 

-(e-m-./?.)^ 1 =f^^^^)'''[(e-^-(''o-^^))sAl^ , 
L ycoshPrJ J L -I/3t 

as depicted in Fig. 3. In particular, if we set Po = Pr, we can rederive one of the equalities 
given in our previous study^^^ as 

_ /cosh^\^^ 
/9o \coshPoJ J 



-y(/3o->/3T) 



)SA 



(54) 
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Fig. 3. Two related nonequilibrium processes of simulated annealing in eq. (53). The same abbrevi- 
ations and symbols are used as in Fig. 1. 



4-2 Quantum annealing and simulated annealing 

Collecting the above results, we can derive several relations between simulated annealing 
and quantum annealing. We combine eq. (31) with eq. (51) for the case starting from /3o = 
and ending at /3t = /3i with the quenched randomness satisfying /J^ = /32 to obtain 



(e '^^ )qa 



1 



N 



, ^^ • (55) 

This equation relates two processes ending at the black dot in the left-lower part of Fig. 4. 
Furthermore, the combination of eqs. (33), (53) and (55) reveals a close relation between four 
annealing processes, classical and quantum, drawn in arrows in Fig. 4. 

In addition, a comparison of eq. (41) with eq. (51) for (3 = /3q, f3p = /3t and /3r = shows 



1 



(e-/5o^)QA 



(cosh /JqF) 



N 



-y(/3o->/3T) 



)SA 



(56) 



jSrj, (cosh/3r J)^^ 

Figure 5 illustrates this relation between the inverse statistics of the performed work in quan- 
tum annealing and the exponentiated pseudo work in simulated annealing. 

The identities derived above constitute a class of new exact relations between classical 
and quantum systems. 

5. Summary 

We have established a number of exact identities and inequalities for nonequilibrium (dy- 
namical) properties of spin glasses in external fields. Essential ingredients in the analysis were 
the gauge symmetry of the Hamiltonian and the Jarzynski equality (or the fluctuation the- 
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Fig. 4. Several relations between simulated annealing and quantum annealing by eqs. (33), (53), and 
(55). The same abbreviations and symbols are used as in Fig. 1. 




Fig. 5. The relation between simulated annealing and quantum annealing in eq. (56) is shown in two 
arrows. 



orem). The results are to be contrasted with our previous study where we discussed spin 
glass systems in the absence of external field and derived identities involving pseudo work for 
the change of the temperature. The present paper deals with genuine work realized by the 
change of the strengths of external fields. The identities and inequalities derived here are rare 
examples of exact results on nonequilibrium properties of spin glass systems and are expected 
to lay a foundation of further analyses of dynamical properties of spin glasses and related 
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problems. 

The Jarzynski equality in the case of transverse field would be a step toward non-adiabatic 
realization of quantum annealing. We remark that it is necessary to repeat the nonequilib- 
rium processes exponentially many times in general to correctly evaluate the nonequilibrium 
average in eq. (47). In this sense we have not solved the problem of exponentially long compu- 
tation time for hard optimization problems in the usual quantum adiabatic evolution but have 
replaced it with exponentially many trials of quick (non-adiabatic) processes. We nevertheless 
expect that the present new point of view would provide a perspective different from previous 
studies, which may help us approach the problem by different methods, eventually leading to 
unexpected results. 

Also, the non-trivial relations between the exponentiated work in quantum and simulated 
annealing may become a valuable tool to compare performance of these two generic algorithms 
for optimization problems. 
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